AD-A191  977  CURRENTS  INDUCED  IN  A  HUNAN  BEING  FOR  ELECTROHAOMETIC 
FIELDS  IB  KHZ-50  HHZ(U)  UTAH  UNIV  RESEARCH  INST  SALT 
LAKE  CITV  0  P  GANDHI  15  JAN  88  ETL-888B2-TR 
UNCLASSIFIED  NBM14-8S-K-8184  F/O  f/it 


one  file  copy 


ETL-88002-TR 


University  of  Utah  Research  Institute 
391  Chipeta  Way,  Suite  C 
Salt  Lake  City,  Utah  84108-1295 


Naval  Medical  R  &  D  Command 
NMC  -  NCR 

Bethesada,  Maryland,  20814-5044 
Attn:  Capt.  W.  M.  Houk 


Final  Report  on  Project 
N  00014 -86 -K  0104 


dtic 

ELECTEI 

UAR07B88 

1 


CURRENTS  INDUCED  IN  A  HUMAN  BEING  FOR 
ELECTROMAGNETIC  FIELDS  10  kHz- 50  MHz 


Principal  Investigator 


CTYv\  Y  Z 

Om~  P .  Gandhi- 


January  15,  1988 


88  2 


0  9 


1  ? 


I  4  **  ||h|k  J 


CURRENTS  INDUCED  IN  A  HUMAN  BEING  FOR  ELECTROMAGNETIC 
FIELDS  lOKHz  -  50  MHz 


inupouraan  '  1 

/ 

^  We  hav^previously  shown  [1 ,  2]  that  vertically  polarized  incident  plane  waves  are  capable  of 
inducing  fairl/  significant  RF  currents  in  a  free-standing  human  being.  Foot  currents  were  found  to  be 
proportional  to  the  frequency  of  incident  radiation  for  the  frequency  band  0-4C^Hz  with  values  as  high  as 
12.7  f$A/(V/m)  measured  for  adult  human  volunteers  at  40  MHz.  This  work  has  pointed  to  the  possibility  of 
large  induced  currents  estimated  to  be  on  the  order  of  627.4  mA  (or  the  ANSI  recommended  [3]  3-30  MHz 
RF  safety  level  -  RMS  incident  electric  field  E  *  1897/Imhz.  with  current  values  peaking  at  780  mA  for  the 
ANSI  recommended  E  «  61 .4  V/m  (1  mW/cm2)  at  40  MHz.  Recognizing  that  the  induced  current  is 
divided  equally  between  the  two  legs  on  its  way  to  the  ground  underneath,  fairly  high  current  densities 
result  in  the  various  cross  sections  of  the  leg  with  concommrtant  high  rates  of  energy  deposition  (SAR).  In 
particular,  due  to  the  predominantly  bony  nature  of  the  ankle  cross  section  forcing  the  RF  current  to  flow  in 
an  effective  9.5  cm2  cross  section  of  the  high  conductivity  tissues,  very  high  current  densities  (J)  and  the 
resulting  SARs  are  set  up.  -Fpr  the  ANSI-recommended  electric  fields,  current  densities  and  SARs 
calculated  for  the  ankle  section  are:  1>for  3-30  MHz  band,  J  -  33  mA/cm2,  SAR  •  182  W/kg  and  2.  for  f « 
40  MHz,  J  ■  41  mA/cm2,  SAR  *  243  W/kg  These  SARs  are  almost  two  orders  of  magnitude  larger  than 

J- 

the  metabolic  rates  of  the  tissues  and  considerably  in  excess  of  the  ANSI  guideline  [3]  of  8  W/kg  for  any 
1  g  of  tissue. 

. ^  The  two  tasks  undertaken  during  the  period  of  this  project  are: 

1 .  Development  of  a  high  resolution  thermal  model  of  the  human  leg  and  its  use  to  estimate 
the  temperature  distribution  as  a  result  of  the  induced  high  current  densities  at  radio 
frequencies. 

2.  Development  of  a  high-resolution,  anatomically-realistic,  inhomogeneous  model  of  man 


and  its  use  to  calculate  the  SAR  and  induced  current  distributions  for  frequencies  to  100  fj 
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TEMPERATURE  DISTRIBUTIONS  IN  7VIE  HUMAN  LEG  FOR 
VLF-VHF  EXPOSURES  AT  THE  ANSI-RECOMMENDED  SAFETY  LEVELS 

This  work  is  detailed  in  Appendix  A  which  has  been  submitted  for  possible  publication  to  IEEE 
Transactions  on  Biomedical  Engineering. 

Using  a  block  model  of  1,532  cubical  cells,  temperature  distributions  are  calculated  for  the  lowest 
21  cm  of  the  human  leg  for  electric  fields  recommended  in  the  ANSI  RF  safety  guideline.  The  thermal 
model  uses  inhomogeneous  volume-averaged  tissue  properties:  blood-flow  rate,  metabolism,  thermal 
conductivity,  specific  heat,  etc.  The  SARs  are  obtained  using  the  impedance  method.  A  modified 
method  of  finite  difference  technique  is  used  to  solve  the  3-D  heat  conduction  equation  for  the  thermal 
model.  Numerical  results  are  obtained  for  RF  currents  at  3  and  40  MHz  projected  for  the  E-fields 
recommended  by  the  ANSI  Standard  (614  and  61 .4  V/m,  respectively)  and  also  for  power  densities  one- 
tenth  of  that  level.  Temperatures  as  high  as  41.6®  are  obtained  for  some  internal  cells  for  the  higher  E- 
fields  while  relatively  moderate  temperatures  on  the  order  of  37°  C  are  obtained  for  the  lower  E-fields. 
Some  of  the  calculated  results  for  the  surface  temperature  have  been  compared  and  found  to  be  in  good 
agreement  with  the  experimental  data  for  initial  rates  of  heating. 

AN  ANATOMICALLY-BASED  MODEL  OF  MAN  AND  ITS  USE  FOR 
SAR  AND  INDUCED  CURRENT  DISTRIBUTIONS 

Similar  to  references  (4]  and  [5],  the  inhomogeneous  model  of  the  human  body  is  based  on  the 
anatomical  cross  sections  given  in  the  book,  ’A  Cross-Section  Anatomy,"  by  Eycleshymer  and 
Schoemaker  (6).  This  book  contains  cross-sectional  diagrams  of  the  human  body  which  were  obtained  by 
making  cross-sectional  cuts  at  spacings  of  about  one  inch  in  human  cadavers.  The  process  for  creating 
the  data  base  of  the  man  model  was  the  following:  first  of  all  a  quarter-inch  grid  was  taken  for  each  single 
cross-sectional  diagram  and  each  cell  on  the  grid  was  assigned  a  number  corresponding  to  one  of  the  14 
tissue  types  given  in  Table  1.  Thus  the  data  associated  with  a  particular  layer  consisted  of  three  numbers 
for  each  square  cell:  x  and  y  positions  relative  to  same  anatomical  reference  point  in  this  layer,  usually  the 
center  of  the  spinal  cored,  and  an  integer  indicating  which  tissue  that  cell  contained.  Since  the  cross- 
sectional  diagrams  available  in  reference  [6]  are  for  somewhat  variable  separations  typically  2.3-2.7  cm,  a 
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new  set  of  equispaced  layers  were  defined  at  1/4"  intervals  by  interpolating  the  data  onto  these  layers. 
Since  the  cell  size  of  quarter-inch  is  too  small  for  ths  memory  space  of  present  day  computers,  the 
proportion  of  each  tissue  type  was  calculated  next  for  somewhat  larger  cells  of  size  one  inch  by  combining 
the  data  for  4  x  4  x  4  *  64  cells  of  the  smaller  dimension.  Without  changes  in  the  anatomy,  this  process 
allows  some  variability  in  the  height  and  weight  of  the  body.  We  have  taken  the  final  cell  size  of  2.62  cm 
(rather  than  1"  or  2.54  cm)  to  obtain  the  whole-body  weight  of  69.6  kg  for  the  model. 

The  finite-difference  tissue-domain  method  described  in  our  earlier  publications  [4,5]  has  been 

used  to  calculate  the  internal  electric  fields  (E)  for  the  model  of  man  subjected  to  vertically-polarized  plane 

waves  at  frequencies  between  20  and  100  MHz.  The  model  is  assumed  to  be  nestled  in  a  rectangular 

volume  of  38  x  26  x  84  or  82,992  cells  of  dimensions  2.62  cm  on  each  of  the  sides.  The  electromagnetic 

fields  are  calculated  using  CRAY  XMP  computer  requiring  cpu  times  on  the  order  of  2.5  -  5  minutes, 

depending  upon  the  frequency.  Both  isolated  and  grounded  conditions  are  used  for  the  calculations. 

— >  — > 

The  local  SARs  are  calculated  from  the  relationship  SAR  «  o  E  •  E  /2p'  where  a  is  the  electrical 
conductivity  and  p  the  mass  density  obtained  for  the  individual  cubical  cell  from  the  volume-averaged 
properties  of  the  respective  tissues  contained  therein.  The  layer-averaged  values  of  SARs  calculated  for 
incident  plane  waves  of  power  density  1  mW/cm2  (incident  E-field  of  61 .4  V/m)  are  shown  in  Figs.  1-10  for 
isolated  and  grounded  conditions  for  frequencies  20-100  MHz.  Lower  frequencies  have  not  been  used 
for  calculations  since  the  SARs  are  considerably  lower  at  these  frequencies. 

The  whole-body-averaged  SARs  for  grounded  and  isolated  man  models  are  shown  in  Fig.  1 1  for 
various  frequencies.  Also  shown  for  comparison  is  the  experimental  data  by  Guy  and  Chou  [7].  Whereas 
these  authors  have  used  scaled,  homogeneous  models  with  conductivity  corresponding  to  2/3  that  for 
muscle,  an  anatomically-based  inhomogeneous  model  has  been  used  for  the  present  calculations.  It  is 
interesting  to  note  that  the  SARs  given  by  the  inhomogeneous  model  are  significantly  higher  than  those 
obtained  for  the  homogeneous  model.  As  will  be  seen  later  (Figs.  13  -  22)  substantial  radio  frequency 
currents  are  induced  in  the  legs.  Because  of  the  large  bone  content  of  the  legs,  an  inhomogeneous 
model  properly  accounts  for  the  lower  effective  conductivities  for  this  region,  which  therefore  results  in 
higher  SARs. 
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From  Figs.  1-10  it  may  be  seen  that  the  highest  SARs  are  calculated  tor  the  sections  through  the 

neck,  the  knee  and  the  ankle  -  the  former  being  due  to  a  smaller  cross  section  and  the  latter  two  due  to  the 

high  bone  content  as  well.  The  layer-averaged  SARs  (or  ankle- ,  knee-,  and  neck  sections  calculated  tor 

grounded  and  isolated  conditions  are  shown  in  Fig.  12. 

We  have  used  the  internal  E-fields  to  calculate  the  local  current  densities  from  the  relationship  J  > 
— ► 

(a + jcjc)  E .  The  z-directed  currents  for  any  of  the  layers  are  obtained  by  summing  up  the  terms  due  to  the 
individual  cells  in  a  given  layer  as  follows: 

1 «  A2  X(°.  (1) 

where  A2  is  the  cross  sectional  area  (2.62  x  2.62  cm2  )  for  each  of  the  cells. 

The  layer-averaged  induced  current  distributions  are  shown  in  Figs.  13-22  for  isolated  and 
grounded  conditions  for  frequencies  of  27-90  MHz.  Since  E-fields  of  61.4  V/m  have  been  recommended 
as  safe  in  the  ANSI  C95.1  -  1982  RF  safety  guideline  J3)  for  the  frequency  band  30-300  MHz,  substantial 
RF  induced  currents  are  obviously  implied  from  the  data  given  in  Figs.  13-22.  A  comparison  is  made  in 
Fig.  23  of  the  foot  currents  calculated  for  a  grounded  model  to  those  that  were  measured  by  Gandhi,  et  al. 
[2]  using  human  subjects.  Somewhat  lower  currents  are  obtained  from  the  model  calculations  than  were 
measrued  experimentally  [2]  This  may  be  due  to  the  relative  crudeness  of  the  present  model.  Finer 
models  using  smaller  cell  sizes  would  provide  better  representation  of  reality  and  may  result 
in  giving  a  better  agreement  with  the  measured  data. 
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Fig.  1 .  Layer-averaged  SAR  distributions  (or  a  grounded  and 
isolated  man  model  at  20  MHz.  All  SARs  are  the  result  of 

1  incident  plane-wave.  The  SARs  (or  the  isolated 

man  model  have  been  multiplied  by  10. 
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Fig  2.  Layer-averaged  SAR  distributions  for  a  grounded  and 
isolated  man  model  at  27  MHz.  All  SARs  are  the  result  of 

inci^ent  P'ane-wave.  The  SARs  for  the  isolated 

condition  have  been  multiplied  by  10. 
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Fig.  3.  Layer-averaged  SAR  distributions  tor  a  grounded  and 
isolated  man  model  at  40  MHz.  All  SARS  are  the  result  of 
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t - y  incident  plane-wave. 
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Fig.  4.  Layer-averaged  SAR  distributions  tor  a  grounded  and 
isolated  man  model  at  45  MHz.  All  SARs  are  the  result  of 
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Fig.  5.  Layer-averaged  SAR  distributions  tor  a  grounded  and 
isolated  man  mode!  at  50  MH*.  All  SARs  are  the  result  of 
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Fig  6.  Layer-averaged  SAR  distributions  for  a  grounded  and 
isolated  man  model  at  55  MHz.  All  SARs  are  the  result  of 

_  rrhN  ...... 

1 - y  incident  plane-wave. 
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Fig.  9.  Layer-averaged  SAR  distributions  lor  a  grounded  and 
isolated  man  model  at  70  MHz.  All  SARs  are  the  result  of 
.  mW 

1 - w-  incident  plane-wave. 
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Fig.  10.  Layer-averaged  SAR  distributions  for  a  grounded  and 
isolated  man  model  at  100  MHz.  All  SARs  are  the 
rrfW 

result  of  1 - y  incident  plane-wave. 
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Fig.  11.  Whole-body  averaged  SARs  lor  a  grounded  and 
isolated  man  model.  Also  shown  for  comparison  is  the 
homogeneous,  isolated  model  data  of  Guy  and  Chou 
17]. 
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Fig.  12.  Layer-averaged  SARs  for  ankle-,  knee-  and  neck 
sections  under  grounded  and  isolated  conditions.  All 

SARs  are  the  result  of  1  incident  plane-wave. 
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Fig  16  Induced  RF  current  distributions  tor  a  grounded  and 
isolated  man  model  at  45  MHz  under  1  V/m  plane-wave 
exposure  condition. 
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Induced  RF  current  distributions  for  a  g 
isolated  man  model  at  48  MHz  under  1  V/r 
exposure  condition. 
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Fig.  19.  Induced  RF  current  distributions  for  a  grounded  and 
isolated  man  model  at  58  MHz  under  1  V/m  plane-wave 
exposure  condition. 
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Fig.  20.  Induced  RF  current  distributions  for  a  grounded  and 
isolated  man  model  at  70  MHz  under  V/m  plane-wave 
exposure  condition. 
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Fig  23.  Comparison  of  RF  foot  currents  calculated  for  a 
grounded  man  model  with  those  measured  for  human 
subjects  [2J .  Ejnc  *  1  V/m. 
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Table  1.  (continued) 

Tissue  properties  used  for  the  man  model. 
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ABSTRACT 


Using  a  block  model  of  1532  cubical  cells,  temperature 
distributions  are  calculated  for  the  lowest  21  cm  of  the  human  leg  for 
electric  fields  recommended  in  the  ANSI  RF  safety  guideline.  The 
thermal  model  uses  Inhomogeneous  volume-averaged  tissue  properties: 
blood-flow  rate,  metabolism,  thermal  conductivity,  specific  heat,  etc. 
The  SARs  are  obtained  using  the  impedance  method.  A  modified  method  of 
finite  difference  technique  is  used  to  solve  the  3-D  heat  conduction 
equation  for  the  thermal  model.  Numerical  results  are  obtained  for  RF 
currents  at  3  and  40  MHz  projected  for  the  E-fields  recommended  by  the 
ANSI  Standard  (614  and  61.4  V/m,  respectively)  and  also  for  power 
densities  one-tenth  of  that  level.  Temperatures  as  high  as  41.6°  are 
obtained  for  some  internal  cells  for  the  higher  E-fields  while 
relatively  moderate  temperatures  on  the  order  of  37°C  are  obtained  for 
the  lower  E-fields.  Some  of  the  calculated  results  for  the  surface 
temperature  have  been  compared  and  found  to  be  in  good  agreement  with 
the  experimental  data  for  Initial  rates  of  heating. 
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1.  INTRODUCTION 

The  present  radio  frequency  RF  radiation  standards  sugg  sted  by 
American  National  Standards  Institute  (ANSI)  and  the  American  Conference 
of  Government  Industrial  Hygienists  are  based  on  the  knowledge  of  whole- 
body-average  specific  absorption  rate  (SAR)  11,2).  For  the  electric 
fields  recommended  in  the  ANSI  C95. 1-1982  RF  Safety  Guide,  fairly  high 
SARs  have  been  projected  for  the  high-water-content  tissues  of  the  ankle 
section.  The  SARs  are  estimated  to  be  on  the  order  of  182  W/kg  for  the 
frequency  band  3-30  MHz  rising  to  a  value  of  243  W/kg  at  40  MHz  for 
E»61 .4  V/m.  The  corresponding  current  through  each  of  the  legs  la 
estimated  to  be  314  mA  and  390  mA,  respectively  13,4).  Our  present 
study  determines  the  temperature  distributions  of  the  lowest  portion  of 
the  leg  under  these  conditions.  For  this  purpose  inhomogeneous 
electromagnetic  and  thermal  models  are  developed  for  the  lowest  21  cm  of 
the  human  leg. 

Temperature  distribution  in  different  parts  of  the  human  body  under 
EV  exposure  is  of  basic  interest  for  radio  frequency  hazards  analyses. 
The  thermal  models  in  existence  to  date  are  capable  of  giving  only 
averaged  temperature  over  large  regions  of  the  body.  The  most  useful 
and  detailed  thermoregulatory  models  used  to  date  are  based  on  a 
cylindrical  thermal  model  described  hy  Stolvijk  (5-7J.  Later  several 
modified  versions  of  this  model  were  developed  by  Wissler  (8,9). 
Lastly,  a  476-block  thermal  model  was  developed  hy  Chatterjee  and  Gandhi 
(10).  This  thermal  model  used  Hagmann's  (11)  180-cell  block  model  of 
man  for  electromagnetic  energy  deposition.  Recently,  a  41,000  cell 
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model  of  wan  has  been  developed  for  SAR  calculations  by  Sullivan,  Gandhi 
and  Taflove  using  finite-difference  tine-domain  technique  112).  In  this 
study  we  have  used  an  Impedance  method  for  detailed  SAR  calculations. 
The  SARs  thus  determined  are  used  In  the  thermal  block  model  to  find  the 
internal  temperature  distribution.  For  temperature  calculations,  the 
lowest  21  cm  section  of  the  leg,  as  shown  in  the  insertion  in  Fig.  2,  Is 
modeled  by  1532  cubical  cells  each  of  size  Ixlxl  cm. 

II.  THE  ELECTROHACNETIC  MODEL  FOR  ENERGY 
DEPOSITION  IN  THE  LEG  SECTION. 

The  model  for  calculations  of  specific  absorption  rates  (SARs)  Is 
developed  with  21  transverse  sections  of  1  cm  thickness  (see  the  Insert 
In  Fig.  2).  For  each  section  the  tissue  types  are  determined  with  a 
resolution  of  1  mm  .  The  individual  sections  of  the  model  are  collected 
from  the  literature  on  the  human  anatomy  (13).  The  model  i6  based  on 
equipotential  concept.  Each  section  is  assumed  to  be  an  equlpotential ; 
i.e.,  no  lateral  flow  of  current  is  supposed  to  exist.  This  concept  can 
be  justified  from  the  fact  that  for  highly-conducting  tissues  (at  VLF- 
VHF  freauencies)  the  potential  gradient  exists  mainly  along  the  axial 
direction  of  the  model.  This  results  in  a  current  primarily  along  the 
axial  direction.  Also,  because  of  the  large  wavelength  of  the  radiation 
no  appreciable  gradient  along  the  transverse  direction  of  the  model 
exists.  As  a  result  the  equivalent  impedance  for  each  section  is 
composed  of  parallel  combination  of  the  impedances  obtained  from  the  ten 
tissue  types  as  6hown  [14,15,16]  in  Table  1. 

Finally,  the  complete  impedance  of  the  model  is  the  series 
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Table  1.  Dielectric  properties  of  the 
different  tissue  types. 


Tissue  types 

Tissue 

type 

No. 

1 

f  - 

Dielectric 

constant 

Ei 

3  MHz 

Conductivity 
o  S/m 

Fat /yellow  rarrov 

1 

33.0 

0.051 

Muscle 

2 

107.0 

0.8 

Red  marrow 

3 

231.0 

0.3 

Tendon /connec t lve 
tissue 

4 

6.5 

o.ooni 

Tendon  and  fat 

5 

20.0 

0.0195 

Bone 

6 

26.0 

0.015 

Skin 

7 

108.0 

0.62 

Blood 

8 

80.0 

1.1 

Synovial  fluid 

9 

95.0 

1.1 

Perlostlum 

10 

95.0 

1.1 

f  -  40  MHz 


constant 

e. 


29.0 

96.0 

138.2 

6.5 

12.11 

15.47 

108.0 

80.0 

95.0 

95.0 


o  S/m 


0.0001 

0.0259 

0.027 

0.62 

1.1 

1.1 

1.1 
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summation  of  the  individual  Impedances  for  each  of  the  sections.  The 
electromagnetic  energy  deposited  in  different  tissues  of  each  section 
can  be  obtained  from  the  corresponding  impedances  and  currents. 
Mathematically  this  may  be  done  as  follows, 


(Oj  "  3“^)  L/[(°i  +  “2  ^A  ] 


1.1 


ij 


where 


i.l 


ii 


fiJ  *  K  +  >  Ei^  Aij/L 


Irpedance  of  the  ith  tissue  in  the  Jth  section 
Admittance  of  the  ith  tissue  in  the  Jth  section 
Area  of  the  ith  tissue  in  the  Jth  section 
Conductivity  of  the  ith  tissue 
Electrical  permittivity  of  the  ith  tissue 
Thickness  for  each  of  the  sections  (assumed  to  be 
in2  m). 


(1) 


(2) 


Power  Pj^  absorbed  in  Z  ■  1^  1^  Re  (Z^) 


*T  TT  °1  Ai1 
10  10 

I  r-  I  4-1- 

i-1  4iJ  1«1  Z 


in 
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(3) 


For  the  complete  model  the  total  values  can  al60  be  written: 
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where. 


■  Total  current  flowing  through  the  leg  section, 

*  Current  flowing  through  the  ith  tissue  of  the  jth  section, 
*=  EM  energy  deposited  in  the  ith  tissue  of  the  jth  section, 

*  Total  impedance  for  the  jth  section, 

*  Total  impedance  of  the  model, 

«  Total  EM  energy  deposited  in  the  model. 


Total  impedance  of  the  model  is  calculated  and  found  to  be 
104  fi  and  85  fi  at  3  MHz  and  40  MHz,  respectively.  These  results  are 
consistent  with  the  experimental  values  (17].  The  6mall  discrepancies 
in  the  impedances  are  because  of  the  discretization  of  the  leg  section 
into  finite  segments. 

The  electromagnetic  energy  deposited  is  thus  calculated  in  the 
individual  sections  for  each  of  the  tissue  types.  The  layer-average 
SAP.s  calculated  for  the  various  sections  are  shown  in  Fig.  1.  The 
figure  indicates  that  the  highest  EM  energy  is  deposited  near  the  ankle 
region  where  highly-conducting  tissues  like  perio6tium,  synovial  fluid 
and  red  marrow  are  present.  It  may  be  noted  here  that  Fig.  1  indicates 
the  average  SAR  values  in  Individual  layers,  while  in  reality  very  high 
SARs  (of  the  order  of  150  W/kg)  in  the  conducting  tissues  in  some  of  the 
sections  (near  the  ankle  region)  are  obtained  with  a  resolution  of 
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Section-average  electromagnetic  energy  deposit 
1  cm  sections  of  the  leg.  I  -  314  mA  at  3  MHz 


5 

1  an  .  Finally,  the  EM  energy  deposition  In  the  block  thermal  model  is 
determined  from  the  energy  absorbed  in  ten  different  tissues  in  each 
section  of  the  electromagnetic  model.  Then  the  heat  conduction  equation 
is  applied  to  find  the  temperature  distribution. 

III.  BLOCK  MODEL  AND  THE  HEAT  CONDUCTION  EQUATION 

Even  though  cylinder  models  have  been  used  in  the  past  (6],  we  have 
formulated  a  block  model  for  the  heat  conduction  equation  in  order  to 
incorporate  the  detailed  SAR  distribution  under  electromagnetic  exposure 
conditions.  Due  to  the  fact  that  localized  electromagnetic  energy 
deposition  is  very  difficult  to  obtain  in  the  cylindrical  model,  the 
block  thermal  model  is  used  to  account  for  such  deposition.  The 
complete  block  thermal  model  is  composed  of  1532  cubical  cells.  The 
weight  of  the  model  is  1.22  kg.  The  cells  are  of  dimension  1  cm  cube 
for  which  the  parameters  are  obtained  from  the  inhomogeneous  volume- 
averaged  tissue  properties  like  thermal  conductivity,  specific  heat, 
blood-flow  rate  and  metabolic  heat  production.  These  parameters  are 
obtained  from  the  physiological  literature  [5,  18-21]  for  13  different 
tissues  as  6hown  in  Table  2. 

The  percentage  of  the  different  tissue  types  in  each  cell  of  the 
model  is  determined  from  the  anatomy  handbooks  [13,22].  The  volume- 
weighted  average  values  of  the  physiological  parameters  in  one 
centimeter  sections  are  determined  by  linear  interpolation.  Total  basal 
metabolic  heat  production  for  the  leg  model  is  found  to  be  0.35  W  and 
the  total  blood  flow  in  the  model  is  determined  to  be  1.988  i/hr.  These 
values  are  consistent  with  the  literature  [6]. 
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Table  2.  Different  tissue  properties  for 
the  thermal  model. 


Tissue  tvpe 

1 

Density 
t’nits  of  103 
kg/m3 

Specific 

heat 

Whr/ °C  kg 

Thermal 
conductivity 
W/°C  m 

Blood 

flow 

rate 

t/min.kp 

Metabolic 

rate 

W/kg 

0.28 

0.00925 

0 

0 

Dry  bone  and 
ligament / tendon/ 
connective  tissue 

1.28 

0.54 

0.69 

0.005 

.02 

Fat /yellow  marrow 

0.85 

0.64 

0.21 

0.02 

0.23 

Muscle 

1.05 

1.05 

0.642 

0.0275 

0.67 

Red  marrow/brain 

1.05 

1.03 

0.528 

0.02 

0.32 

Tendon /I i g ament / 
connective  tissue 

1.05 

0.64 

0.21 

0.01 

0.004 

Dry  bone  and 
yellow  marrow 

1.2 

0.54 

0.69 

0.01 

0.02 

Fat  and  skin 

0.925 

0.85 

0.275 

0.063 

0.5 

Tendon/ ligament / 
connective  tissue 
and  muscle 

1.05 

0.85 

0.43 

0.019 

0.34 

Tendon /ligament/ 
connective  tissues 
and  fat 

0.95 

0.64 

0.21 

0.015 

0.015 

The  transient  heat  conduction  equation  {10]  which  is  solved  for  the 
theraal  response  of  the  block  model  is  given  by  (5). 

JT 

P4  Y r*  -  V  (k.V  T.  )  ♦  h  .  +  h  -  h  .  -  h  .  -  h  .  +  Q..  (T.  -  T,  ) 

i  1  at  i  irl  ewi  ci  radi  convi  Hbi  v  b  1' 

(5) 

The  significance  of  the  different  parameters  16  discussed  in  the 
appendix. 

The  steady-state  temperatures  obtained  in  the  peripheral  cells  of 
the  model  are  weiphted  average  temperatures  over  the  volumes  of  the 
cell6.  The  peripheral  cells  contain  a  substantial  amount  of  air  and  as 
a  result  the  volume-average  thermal  conductivities  of  the  cells  are 
found  to  be  very  small  which  prevented  the  normal  heat  flow  from  the 
internal  section  of  the  model.  This  is  overcome  by  assuming  the  thermal 
conductivities  of  the  cell  to  be  averaged  over  the  tissue  types  of  the 
cell.  The  same  approach  is  adopted  for  the  specific  heat  of  the 
peripheral  cells.  The  surface  area  of  the  thermal  block  model  is  larger 
than  the  Dubois  surface  area  of  a  standard  man.  For  this  reason,  some 
multiplier  factors  are  introduced  which  account  for  the  larger  surface 
area  of  the  thermal  model  by  decreasing  the  effective  surface  area  for 
the  heat  exchange  terms  to  values  given  in  [23,24).  For  this  purpose 
the  model  is  considered  in  two  separate  sections,  one  is  the  lower 
section  of  the  leg  (up  to  the  ankle)  while  the  other  is  the  foot.  The 
constants  are  taken  as  the  Tatio  of  the  Dubois  surface  area  to  the  model 
area  for  each  of  these  sections. 
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IV.  SOLUTION  OF  THE  HEAT  CONDUCTION  EQUATION  AND  XESULTS 


The  heat  conduction  eauatlon  is  solved  in  rectangular  coordinates 
hy  a  method  which  is  a  modification  of  the  Douglas  and  Rachford  (25) 
that  achieves  higher-order  accuracy  of  the  Grank-Nicholaon 
formulation.  The  unconditionally  6table  system  similar  to  that  of 
Douglas  and  Rachford  method  had  the  advantage  of  using  larger  time 
steps,  while  the  unknown  temperatures  are  computed  explicitly  at  the 
advanced  time  level,  similar  to  that  in  the  Grank-Nicholson  formulation 
having  higher-order  accuracy.  As  a  result  with  reduced  computational 
time  reasonable  accuracy  is  attained.  The  SOR  method  is  used  for  rapid 
convergency  of  the  steady-state  solution.  The  equation  is  solved  using 
the  HP1000  computer. 

For  6teady-state  temperature  distributions  under  basal  condition, 
sweating  is  neglected.  The  calculated  layer-average  basal  temperature 
distributions  are  shown  in  Fig.  2.  Heat  transfer  at  the  bottom  of  the 
model  is  assumed  to  occur  through  a  shoe  having  high  insulation  (23). 
The  lower  temperature  region  near  the  bottom  of  the  foot  i6  because  of 
the  extremity.  The  gradient  inversion  near  the  layer  4  in  Fig.  2  is 
because  of  the  fact  that  the  presence  of  the  shoe  is  considered  up  to 
the  layer  3.  Also,  this  physical  dimension  of  the  model  near  these 
sections  sharply  reduces;  as  a  result,  there  are  more  peripheral  cells 
containing  a  large  percentage  of  air.  At  the  top  section  of  the  model 
heat  transfer  is  assumed  to  occur  only  through  blood  at  35.8°  C.  The 
sudden  sharp  fall  in  the  upper  section  may  be  due  to  the  truncated 
boundary  of  the  model. 


120  mA 


After  obtaining  basal  temperatures,  temperature  distribution*  under 
electromagnetic  energy  depositions  at  one-tenth  of  the  ANSI-recommended 
power  density  levels  is  determined.  The  current  under  this  condition  is 
100  mA  at  3  KHz.  Layer-average-temperature  rise  with  respect  to  the 
corresponding  basal  values  is  al6o  shown  in  Fig.  2.  Consistent  with 
some  preliminary  experiments  with  a  human  subject,  sweating  is  neglected 
under  these  conditions  as  well.  It  is  observed  that  at  100  mA  current 
(f-3  MHz)  the  maximum  rise  In  temperature  is  found  to  be  0.6°C  near  the 
ankle  region,  while  at  120  mA  current  (f-3  MHz)  the  rise  in  temperature 
is  found  to  be  about  1.0°C.  The  heat  balance  for  a  number  of  RF  current 
flow  conditions  Is  shown  in  Table  3. 

Temperature  distributions  are  also  obtained  for  electromagnetic 
energy  depositions  projected  for  the  present  ANSI  Standard  [1].  The  RF 
currents  projected  for  each  of  the  legs  are  314  mA  at  3  MHz  and  390  mA 
at  40  MHz,  respectively  (3].  For  these  conditions,  sweating  is  included 
for  the  peripheral  cells.  Also,  the  Increase  in  the  blood-flow  rate  and 
basal  metabolism  is  taken  into  account.  The  heat  loss  by  sweating  in 
the  ith  cell  is  represented  by  (6,23,24] 

SV  “  he  Aeff  (PWci  *  PVamb)  (6) 

hp  -  Evaporative  coefficient  of  heat  transfer 
Pw  t  ■  Partial  pressure  of  saturated  water  vapor  on  clothing. 

^Wamb  *  Partial  pressure  of  water  vapor  (at  50X/RH  at  30®  C.) 

Also,  the  increase  in  blood-flow  rate  for  the  internal  tissues  at 
higher  temperatures  (above  39°C)  is  given  by  the  linear  relation  as 


Table  3.  Heat  halanre  table  for  banal  and  KM  energy 
deposition  condition  In  the  leg  section. 


Qdllat  “  Bf  for  T1<  39 

Qdilat  -  ®f  [1+0.8  (Ti  -  39)]  for  39°C  <  Tj  <  44°C 

^dilat  “  5  Bf  for  Ti  >  44#c  (?) 

where 

^dilat  “  Increased  blood  flow  1/ hr 
Bf  ■  Basal  blood  flow  1/hr 

From  physiological  considerations,  upper  limits  have  been  taken  for 
both,  the  heat  loss  by  sweating  and  the  Increase  in  blood-flow  rates. 
The  maximum  rate  of  sweating  depends  on  the  water  vapor  pressure  In  air, 
water  vapor  pressure  at  the  skin  and  the  rate  of  convective  heat 
transfer.  The  blood-flow  rates  in  the  peripheral  cells  are  assumed  to 
saturate  at  somewhat  higher  values  than  the  internal  tissues  and  are 
taken  to  be  ten  times  the  basal  values.  To  find  the  saturation  values 
integrated  output  from  skin,  warm  and  cold  receptors  are  determined  [6]. 

The  heat  conduction  equation  is  solved  after  introducing  the  heat 
loss  by  sweating  and  increase  in  blood-flow  rate  and  metabolism.  Table 
3  indicates  that  a  reasonably  good  heat  balance  occurs  for  all  of  the 
cases  that  are  considered.  The  average  rise  in  temperature  for  314  mA 
(at  3  MHz)  and  390  mA  (at  40  MHz)  in  each  section  of  the  model  is  shown 
in  Fig.  3.  The  figure  indicates  that  for  the  prescribed  currents  the 
highest  average  increase  in  temperature  is  3.8°C  and  4.85"C  for  314  mA 
and  390  mA,  respectively,  near  the  ankle  region.  It  may  be  mentioned 
here  that  from  Fig.  1  the  Di  energy  deposition  is  found  to  be  highest  in 
these  sections  as  well. 


390  mA 


For  clarity  of  presentation,  temperature  distribution  in  one  of  the 
sections  near  the  ankle  region  is  shown  in  Fig.  4(b).  For  comparison, 
the  basal  temperature  in  that  section  is  also  given  in  Fig.  4(a).  The 
hot  spots  are  found  to  exist  in  the  internal  region  of  the  section. 
This  is  because  of  the  presence  of  highly-conducting  tissues  like 
synovial  fluid,  periostium,  red  marrow  (present  inside  small  bones  126)) 
and  muscle  which  results  in  high  EM  energy  in  these  tissues.  Because  of 
the  fact  that  some  of  these  tissues  have  thermal  conductivities  which 
are  not  commensurately  high  (see  Table  2)  and  some  are  placed  in  the 

interior  region  (e.g.,  posterior  tibial  muscle  in  Fig.  4(a,b)),  the 

generated  heat  i6  diffused  to  the  surrounding  cells  at  a  lower  rate. 
Therefore,  this  results  in  higher  temperature  under  steady-state 
condition.  Also,  a  smaller  physical  cross  section  for  thi6  part  of  the 
model  (ankle  region)  Is  partially  responsible  for  high  SAR  values.  It 
may  be  observed  that  the  temperature  at  the  frontal  part  of  the  section 
in  Fig.  4(b)  is  higher  than  that  in  the  rest  of  the  peripheral  area. 
This  is  because  this  portion  of  the  section  is  composed  of  highly- 
conducting  tissues  like  periostium  (surrounding  the  tibia  and  fibula), 
synovial  fluid  (present  between  the  tibia  and  fibula),  anterior  tibial 
muscle  and  red  marrow  (present  inside  the  bones).  As  a  result,  high 
SARs  are  deposited  in  this  part  of  the  section  ’Jhich  culminates  in 

increased  temperature.  The  high  rise  in  temperature  in  these  sections 

causes  the  blood  flow  to  increase  appreciably.  The  increased  blood  flow 
due  to  vasodilation  at  314  mA  (f  -  30  MHz)  is  shown  in  Fig.  5.  It 
indicates  that  blood  flow  Increases  to  about  six  to  6even  times  its 
basal  values  in  these  section. 


\  Posterior 
tibial  muscle 

n  and  fat  layer 

(a)  For  basal  condition  without  RF  current. 


(b)  At  1  '  314  mA  (f  ■  3  MHz) 


Contours  of  constant  temperature  (°C)  for  Section  10 
of  the  model  near  the  ankle  region. 


It  is  observed  from  Fig.  3  that  the  average  temperatures  for  the 
lowest  sections  of  the  model  are  slightly  less  than  the  basal 

temperatures.  This  is  because  of  the  presence  of  sweating  which  results 
in  •  drop  in  the  skin  temperature.  This  decrease  in  the  skin 

temperature  as  a  result  of  sweating  has  been  reported  previously  in  the 

literature  [27,28].  The  curves  shown  in  [27]  indicate  that  under 

sweating  condition  the  mean-skin-temperature  of  the  subject  can  drop  to 
as  low  as  30°C  from  a  nonsweating  temperature  of  34.5°C. 

IV.  TRANSIENT  SOLUTION  AND  COMPARISON  WITH  EXPERIMENTS 

In  order  to  compare  the  behavior  of  the  model  with  that  of  the 
recent  experiments  [29],  transient  response  of  the  model  is  determined 
at  5  and  40  MHz  for  three  minutes  and  eight  minutes  of  Irradiation, 
respectively.  The  point  taken  for  comparison  is  at  the  front  of  the  leg 
at  the  ankle  section  (see  inserts  of  Figs.  6  and  7)  which  corresponds  to 
the  region  of  maximum  superficial  heating  [29]  under  conditions  of  RF 
current  flow.  Since  the  initial  temperature  buildup  did  not  produce  any 
observable  sweating,  the  thermal  model  also  did  not  consider  the 
sweating  during  the  first  few  minutes.  Since  the  experiments  of  eight 
minutes  duration  (limited  in  time  by  physical  discomfort)  indicates  an 
initiation  of  sweating  after  about  seven  minutes  [30],  we  too  have 
introduced  such  a  loss  in  our  model.  A  delay  in  the  sweating  may  be 
because  the  full  response  of  the  sweat  glands  requires  a  finite  time 
interval  to  8ttain  its  full  activation  level.  The  results  of  the 
calculations  under  transient  conditions  are  shown  in  Figs.  6  and  7. 
Since  the  exposure  time  in  Fig.  6  is  less  than  6even  minutes,  sweating 
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is  not  considered  in  the  calculations.  It  may  be  noticed  in  Fig.  7  that 
the  surface  temperature  attains  a  maximum  value  of  about  42.6°C,  then  it 
gradually  reaches  to  it6  steady-state  value.  The  sharp  drop  in 

temperature  after  about  12  minutes  is  due  to  Increased  blood  flow  and 
sweating. 

VI .  MODEL  DRAWBACKS 

It  should  he  realized  that  the  model  has  some  deficiencies  in 

different  aspects.  The  cell  6lze  of  the  thermal  model  16  taken  to  1  cm 

cubes,  which  indicates  the  volume  average  properties  of  the  tissues. 
For  this  reason,  6ome  of  the  tissues  such  as  periostium,  synovial  fluid, 
skin,  which  are  considered  for  the  electromagnetic  model  (with  1  mm 
resolution),  cannot  he  taken  into  account  properly  in  the  thermal 
model.  Again,  lack  of  proper  modeling  of  the  skin  (which  is  responsible 
for  the  total  heat  exchange  of  the  model  with  environment)  has 
introduced  6ome  limitations  on  the  model  performance.  It  is  observed 
that  most  of  the  peripheral  cells  contained  more  than  75  percent  air 
which  resulted  in  lowering  of  the  overall  temperature  of  these  cell6. 

The  effect  on  heat  flow  because  of  this  deficiency,  however,  is 
corrected  by  taking  the  values  of  heat  conductivities  and  specific  heats 
corresponding  only  to  the  volume-averaged  properties  of  the  biological 
tissues  in  the  peripheral  cells.  Since  we  used  1  mm  resolution  in  the 
cro6s-sectional  planes  for  the  electromagnetic  model,  the  detailed  SAR 
calculations  are  relatively  free  from  the  large  cell  size  problems.  As 
it  is  not  possible  to  obtain  the  anatomical  cross  sections  with  less 
than  1  cm  spacing,  the  electromagnetic  model  is  obtained  by  connecting 
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21  equivalent  Impedances  In  aeries  with  each  other.  It  has  been 
mentioned  earlier  that  there  16  a  slight  variation  of  the  Impedance 
obtained  from  this  model  as  compared  to  the  experimental  results 
obtained  by  Guy,  et  al.  [17].  This  Is  because  of  the  finite  thickness 
of  the  various  sections  used  for  the  model. 

VII.  DISCUSSION 

This  paper  describes  an  inhomogeneous  block  thermal  model  of  the 
lowest  21  cm  of  the  human  leg  consisting  of  1532  cells  each  1  cubic  cm 
in  6ize  •  The  object  of  this  model  is  to  study  the  temperature 
distribution  In  the  lowest  part  of  the  leg  for  conditions  of 
electromagnetic  exposure  prescribed  by  the  present  ANSI  Standard.  Two 
analyses  are  performed  in  particular.  These  are  for  314  mA  at  3  MHz, 
and  390  mA  at  40  MHz  —  currents  that  may  be  6et  up  in  the  human  leg 
[3,4]  for  E-flelds  recommended  by  the  present  ANSI  Standard.  Fig.  3 
indicates  that  average  Increases  in  temperature  for  the  ankle  section 
under  these  conditions  are  3.8°C  and  4.85CC,  respectively.  Both 
terperature  profiles  indicate  much  similarity.  The  region  of  high 
temperature  is  near  the  ankle.  The  blood-flow  rate  increases 
substantially  under  these  conditions  (Fig.  5). 

The  heat  supplied  by  the  blood  is  found  to  vary  under  these 
conditions.  At  314  mA  (f-3  MHz)  0.1  Watt  of  heat  is  taken  by  the  blood 
In  contrast  to  1.4  Watt  heat  supplied  by  the  blood  under  basal 
condition.  At  390  mA  it  was  observed  that  3.01  Watt  of  heat  is  supplied 
to  the  blood  from  the  leg  model  under  consideration.  The  increase  in 
metabolism  is  also  indicated  in  Table  3.  Studies  are  also  made  to  find 
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the  temperature  distribution  of  the  model  with  EM  power  deposition  for 
Incident  power  densities  roughly  one  tenth  of  those  given  by  the  present 
ANSI  guideline.  As  observed  by  a  subject  experimentally,  sweating  is 
neglected  under  these  conditions.  The  RF  currents  used  are  100  *A  at 
and  120  mA  at  3  MHz,  respectively.  The  section-average  temperatures 
rose  about  0.6*C  and  1.0*0,  respectively,  for  the  ankle  region.  These 
results  are  in  good  agreement  with  the  experimental  data  (29).  Lower 
temperature  increases  are  calculated  for  the  rest  of  the  leg  sections. 
Transient  solutions  of  the  model  are  compared  with  the  experimental  data 
for  the  initial  Increases  in  temperature  near  the  ankle  region. 
Numerical  results  at  399  mA  indicate  s  sharp  rise  in  temperature  up  to 
42.6°C  after  about  10  minutes  of  irradiation.  The  results  are 
consistent  with  the  experimental  data.  Increased  blood  flow  and 
sweating  causes  the  temperature  to  eventually  attain  a  steady-state 
value  of  39.3°C.  As  the  behavior  of  the  sweat  glandB  under  transient 
conditions  is  not  exactly  known,  the  curve  shown  in  Fig.  7  may  have  to 
be  modified  somewhat,  according  to  the  performance  characteristics  rc 
the  6weat  glands. 
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APPENDIX 

The  Bioheat  Equation 

The  transient  bioheat  equation  [10)  vhi fh  Is  solved  for  the  thermal 
response  of  the  block  model  Is  given  by 


1 

P,e.  -r- —  -  V  (K.  VT. )  +  h  .  +  h  .  -  h  .  -  h  ..  -  h  ,  (A-l) 

i  1  3t  1  1 J  ml  eml  el  radl  convl  ' 


*  \i  (Tb  -  V 


Tj  «  Instantaneous  temperature  of  the  cell  (°C) 

«  Mass  density  of  the  1th  cell  (kg/tr^) 

Cj  ■  Specific  heat  of  the  1th  cell  [Whr/(#C  kg)] 
kj  «  Thermal  conductivity  of  the  cell  [W/(°C  m)] 

2 

hml  *  Metabolic  heat  generation/unit  volume  (W/m  ).  (Varies 
with  the  rise  of  temperature.) 


2 

bemi  *  Electromagnetic  energy  deposltion/unlt  volume  (W/m  ) 

hg^  ■  Evaporative  heat  dissipation  in  the  peripheral 
cells/unit  volume  (Includes  sweating  and 
Insensible  perspiration)  (W/mJ) 

h  ,  h  *  Radiative,  convective  heat  loss  from  the  peripheral 

1  cells/unit  volume  (W/m  ) 

°M  “  Product  °f  mass  ^ow  rate»  density  and  specific 

beat  of  capillary  blood  in  the  tissue  [W/(m3  *C)] 

*  Temperature  of  arterial  blood  entering  the  tissue 
assumed  to  be  at  35.8°C. 
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It  is  assumed  that  heat  exchange  between  blood  and  tissue  occurs  in 
the  capillaries  only  15,10).  The  boundary  conditions  are  introduced  in 


the  peripheral  cells  by  the  loss  terms  as  hrad,  hconVi,  and  hfil.  The 
empirical  equations  for  these  loss  terms  have  been  obtained  from  studies 
by  Fanger  (23)  and  Gagge  and  Nishl  1 24 )  for  clothed  body.  These  are 
given  in  the  following: 


hr«dl  '  *eff  'S  l(Tct+  273 )‘  *  (*,  +  2?3)4]  (K  c«l/hr)  U-2) 

\onvj  •  *D>,  fc,  \  [Tcl  -  Til  <*  (*-3) 

hpD  *  0.35  ADu  (1.92  Tj  -  25.3  -  Pfl)  (K  cal/hr)  (A-4) 

h.l  ‘  V,  <>.1)*T  w/»3  <*-5) 

Tct  ‘  T.  *  Fct  (Ti  -  V  ’C  <*-» 


F  0  '  I„/I„  +  I  ,) 

c£  o  o  cl 


Aef f  *  fef f  fcl  ADu 


C  A— 7 ) 

(A-8) 


where , 

Agjf  *  Effective  area  of  the  clothed  body. 

feff  m  Fatio  of  the  effective  radiation  area  of  the  clothed  body 
to  the  surface  area  of  the  clothed  body  »  0.65  (for 
seated  person). 

f  j  -  Patio  of  the  surface  area  of  the  clothed  body  to  the 
nude  body. 

Aj^  ■  The  Dubois  area. 

h^ib  ■  Basal  metabolic  heat  production. 


c 

p« 

rc£ 

AT 

rib 

Ta 

rcl 

I 


c£ 


Convection  coefficient. 

2,05  (Tct  “  V*25*  <K  c«l/B>2-hr-*C) 

-8  -2  -4 

Stefan-Bol troan  constant  ■  5.67  x  10  Wo  K  . 
Emissivlty  of  skin  •  1. 

Partial  pressure  of  water  vapour  in  the  ambient  air 
Clothing  temperature. 


Temperature  of  the  ith  cell  under  basal  condition. 
Ambient  temperature  (30*C). 

Thermal  efficiency  factor. 

Thermal  resistance  of  operating  environment. 
Effective  insulation  of  clothing. 
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